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Abstract 

The new model of modified F[R) gravity theory with the function 
F{R) = R + (a/ 7 ) arcsin( 7 i?) is suggested and investigated. Constant 
curvature solutions corresponding to the extremum of the effective 
potential are obtained. We consider both the Jordan and Einstein 
frames, and the potential and the mass of the scalar degree of freedom 
are found. It was shown that the de Sitter space-time is unstable but 
the flat space-time is stable. We calculate the slow-roll parameters e, 
rj, and the e-fold number of the model. Critical points of autonomous 
equations for the de Sitter phase and the matter dominated epoch are 
obtained and learned. 


1 Introduction 

To describe the inflation and the present time universe acceleration one may 
try to modify the Einstein-Hilbert (EH) action of general relativity (GR). We 
analyze here the particular case of the F{R) gravity model by replacing the 
Ricci scalar R —)■ F{R) in EH action. It is known that F{R) gravity models 
can describe the evolution of universe without introducing Dark Energy (DE) 
n, 0,0. In such models the cosmic acceleration is due to modihed grav¬ 
ity, and F(i?)-gravity models may be an alternative to A-Cold Dark Matter 
(ACDM) model because of the new gravitational physics. In ACDM model 
the cosmological constant A is introduced but there is a problem with the 
explanation of the smallness of A. 

The hrst successful models of F(i?)-gravity describing inflation were given 
in a, 0 , i, 0 , i, i, [To] . Some F(i?)-gravity models were discussed in 
0 , 0 , n, ca, 1131, iia. eg. iig. ini and in many other publications. It 
should be noted that F(R)-gravity is the phenomenological model that can 
describe the evolution of universe and gives the self-consistent inflation. It 
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should be mentioned that models 0.0.131 represent the possibility to unify 
DE with inflation in a realistic way. 

In this paper we introduce a new model of F(R) gravity that is a model 
for inflation and inflationary parameters are studied. 

In Sec.2, we formulate a model of F{R) gravity with the dimensionless 
parameter a and the parameter 7 with the dimension of (length)^. We hnd 
the constant curvature solutions corresponding to the de Sitter space-time. 
In Sec.3, the potential and the mass of the scalar held are obtained in the 
Einstein (scalar-tensor) form of the model. The graphs of the functions 
0 ( 7 /?), and m‘^{jR) are represented. It is shown that the de Sitter 

phase is unstable and the solution with the zero curvature scalar is stable. 
We calculate the cosmological parameters e, 7 , and the e-fold number of the 
model in Sec.4. In Sec.5 critical points of autonomous equations are studied. 
We estimate the function m(r) characterizing the deviation from the ACDM 
model. The Sec .6 is devoted to a discussion. 

The Minkowski metric 7 ^,^=diag(—1,1,1,1) is used and c=h=l is as¬ 
sumed. 


2 The Model 

We consider the modihed F(i?)-gravity model with the function 

F(R) = R + — arcsin( 7 i?), ( 1 ) 

7 

where a is dimensionless parameter and 7 possesses the dimension of (length)^. 
The action of the model in the Jordan frame is given by 

S = ^ j d*x^F{R), (2) 

were k = Mp/, Mpi is reduced Planck’s mass, and we do not include the 
matter Lagrangian density. At a = 0 one arrives at the EH action. According 
to Eq.(l) F{0) = 0, so that R = 0 corresponds to the flat space-time without 
cosmological constant. As GR passes local tests, we assume that the second 
term in Eq.(l) is small at low curvature regime, jR 1, a < 1. We consider 
the description of the universe evolution in this model. 
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Classical stability takes place if F'{R) > 0 [T] so that we require 


a + Jl- (7i?)2 

F'{R) = / ' > 0. (3) 

\/i - bRr 

Inequality (3) is satisfied for a > 0, yi? < 1. Quantum stability leads to 
inequality F'\R) > 0 |1], which reads 


F''{R) = 


a-f^R 


11 - (7-R)^]’'' 
and Eq.(4) is also satishes at a > 0, yi? < 1. 


> 0 , 


( 4 ) 


2.1 Constant Curvature Solutions 

Equations of motion [TH], in the case when the Ricci scalar Rq is a constant, 
give the equation as follows: 

2F{Ro) = RoF'iRo). (5) 

It should be noted that constant curvature solutions of Eq.(5) correspond to 
the extremum of the effective potential. From Eqs.(l), (5) we hnd 

2a^l - (yRo)^ arcsin(yRo) = 7 -Ro - \/l - (yRo)^^ • (6) 

The solution i?o = 0 to Eq.(6) corresponds to the flat space-time. From 
Eq.(6) we obtain the parameter a as a function of a; = 'yRo: 

xVl - 

® =- o /I-2 -^ 

X — 2\/l — x‘^ arcsm x 

The plot of the function a{x) is given by Fig.l. We obtain numerically 
solutions to Eq.(6) for different parameters a which are given in Table 1. 
Solutions given in Table 1 correspond to the Schwarzschild-de Sitter space- 
time and the solution R = 0 corresponds to the Minkowski space-time. The 
inequality F'(Rq)/F"{Rq) > Rq [19] assures that the constant curvature 
solutions describe DE which is future stable. This requirement leads to the 
equation 

(l — x^'^ (a + Vl — x^^ > ax^. ( 8 ) 
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Figure 1: The function a versus 'jRo corresponding to constant curvature 
solutions. 

The solution i?o = 0 obeys Eq.(8) and, therefore, it is stable. The solutions 
given in Table 1 do not satisfy Eq.(8), and they lead to unstable de Sitter 
space-time. We will show that these solutions describe inflation and corre¬ 
spond to the maximum of the effective potential in the Einstein frame. As a 
result, even without matter, the model of Eq. (1) mimics DE and gives new 
non-trivial potential in the Einstein frame. The general study of de Sitter 
space-time solutions in F{R) gravity models was done in [I^. Thus, we show 
that in our specihc model based on Eq. (1) there are many solutions cor¬ 
responding to the de Sitter space-time depending on the parameter a. This 
allows us to analyze the dependence of different universe scenarios on the 
parameter a. 
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Table 1: Constant cnrvatnre solntions 


a 0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

X 0.997 

0.992 

0.987 

0.982 

0.978 

0.974 

0.971 

0.968 

0.965 


3 The Scalar-Tensor Formulation 

One can investigate the theory in the Einstein frame to get the scalar-tensor 
form of the theory. For this pnrpose we perform the conformal transformation 
of the metric [20] 




F\R)g,. = 


-S'/.: 

\/i - (iRY 


(9) 


Then the Lagrangian density in Einstein’s frame becomes 

R = ^R - - E(0), (10) 

where the Ricci scalar R in Einstein’s frame is evalnated in new metric (9). 
The scalar held 0 is given by 


0(R) 


\/2k 


\nF'{R) 


a + ^1- {'jRf' 


( 11 ) 


and the potential V (0) is 


V{R) 


RF\R) - F{R) 
2k^F'^{R) 


1 ^ 1 - (7^)2 [yR - 


27^2 (^a + 


{'jRy arcsin(7R) 

w)' 


( 12 ) 


The plots of the fnnctions K(j) and versns yi? are given in FigJ2] and 

FigJ21 The plot of fnnction y/t2R versns K<p for the valne a = 0.9 is repre¬ 
sented by FigJU It follows from Eq.(12) that the condition for an extremnm 
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Figure 2: The function ncj) versus 'yR. 


of the potential dV/dR = 0 leads to Eq.(5). Therefore, the constant curva¬ 
ture solutions to Eq. (6) correspond to the extremum of the potential. Thus, 
the potential function (12) has the minimum at i? = 0 and the maximum 
given by Table 1. The flat (Minkowski) space-time {R = 0) is the stable 
state and the states with the curvatures given in Table 1 are unstable. 

From Eq.(12), we obtain the mass squared of a scalar degree of freedom 
(scalaron) 

2 dV _ 1 / 1 R 4F{R) \ 

3 \F''{R) ^ F'{R) ~ F'\R)) 


^/l — x'^ 

r (1 

— -I- a 

4\/r 

- x'^ iyx + a arcsin x)' 

37 

-ax 

aF \/\ — x^ 


a + \/l — x"^ 

2 


The plot of the function versus x = 'yR is represented in Fig.5. One can 
verify that < 0 for the constant curvature solutions given in Table 1, and. 
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Figure 3: The function versus 'yR. 


therefore, these solutions correspond to unstable states. The stability of the 
de Sitter solution in F(R) gravity models was hrst studied in [19]. For small 
value of the parameter 7 the mass m^j, is big and, as a result, corrections to 
Newton’s law are negligible. 

Corrections of F{R) gravity model will be small comparing to GR for 
R^ Rq {Rq is a curvature at the present time) if the relations hold [T], 

I F{R) - R |< R, I R'(R) - 1 |< 1, I RR"(R) |< 1. (14) 

From Eqs.(3), (4) and (14) we obtain 

I aarcsina: [-C x, a <C Vl — x^, ax^ <C (^1 — x^^ ^ , (15) 

where x = 7 R. Eqs.(15) are satisfied if 0 < a < 0.9, x < 0.43. 
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Figure 4: The function versus ncj) {a = 0.9). 

4 The Slow-Roll Cosmological Parameters 

Consider the slow-roll parameters |21] 



( 16 ) 


The slow-roll approximation occurs if conditions \ri{(j))\ <C 1, e(0) -C 1 hold. 
When e(0) < 1 the acceleration of universe takes place and inflation ends if 
e, \ri\ ^ 1. At |r7(0)| < 1, e(0) < 1 the kinetic energy of scalar field is small 
during inflation. We obtain the slow-roll parameters from Eqs.(ll)-(13),(16) 

1 \RF'{R)-2F{R)Y 
^ “ 3 [ RF'{R) - F{R) _ 
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Figure 5: The function '^m'^ versus 7 -R. 


1 

” 3 
r] = 


X 1 

(vi - 

x^ — aj 

1 + 2a\/l — x^ arcsin X 

a 1 

[yi — x"^ arcsin x — xj 

1 


F'\R) + F"{R) [RF'{R) - AF{R)] 
F''{R) [RF'{R) - F{R)] 


+ (1 — + a\/l — (2 — 5x^ — Aax arcsin x) 


?txa? (^x — \/l — arcsin x^ 


(17) 


(18) 


where x = 'yR. The plots of the functions e, rj are represented in Figs.6 and 7. 

The equation e = 1 has the solution x ~ 0.876 for a = 0.9 and x ~ 0.892 for 
a = 0.7. It follows from Fig.6 that at 1 > 'jR > 0.876 (a = 0.9) the inequality 
e < 1 holds. The equation |? 7 | = 1 is satished at x ~ 0.633, x ~ 0.653 and 
X ~ 0.971 (see Fig.7). At 0.653 > 'jR > 0.633 and at 1 > 'jR > 0.971 
(a = 0.9), we have the result l^l < 1. Therefore, the slow-roll approximation. 
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Figure 6 : The function e versus 7 -R (a=0.9). 


e < 1 and |? 7 | < 1, takes place at 1 > 7 -R > 0.971 (a = 0.9). The calculation 
of inflationary parameters allowed us to hnd the range of curvature when the 
slow-roll approximation of inflation takes place. 

We can obtain the age of the inflation by calculating the e-fold number 


PH 






(19) 


From Eqs.(ll),(12) we find the number of e-foldings 




3a^ 


r^o 


X (^x — \/l — arcsin x^ dx 


(1 — x^) (^a -|- \/l — x'^'^ \/l — {x -I- 2a arcsin x) — 


ax 


( 20 ) 

were Xend = iRend- It should be noted that integral ( 20 ) is singular at x = 
7 R = 1 and, therefore, the value N^, is sensitive for the range of integration. 


10 












Figure 7: The function t] versus 7 /? (a=0.9). 


We obtain the necessary amount of inflation W ~ 40 at Xend = 0.999 and 
xq = 0.95 for a = 0.9. Thus, the model under consideration describes the 
inflation and can give the reasonable age of the inflation. The tensor-to- 
scalar ratio is given by the relation [21] r = 16e. The PLANCK experiment 
gives the restrictions ra r < 0.11. From this bound we obtain from Eq. 
(17) the range of curvatures 0.95769 < yi? < 0.97223 (a = 0.9) realizing 
the tensor-to-scalar ratio restriction. We note that r = 0 at yi? ~ 0.965 
(a = 0.9). 
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5 Critical Points of Autonomous Equations 


To investigate critical points of equations of motion, it is useful to introduce 
the dimensionless parameters [23] which become 


F'{R) axx 

HF'{R) H + aj [1 — x"^) ’ 


F{R) {x + a arcsin x) \/l — x^ 

6F'{R)H'^ + \/l — 



2 , 


_ RF''{R) _ __ 

(1 — a;^) (a + y/l — x‘^'^ ’ 

RF'{R) _X 3 _ x(a + Vl - a;2) 
F{R) X2 vT — x‘^ {x + a arcsin x) ’ 


( 21 ) 


( 22 ) 


where x = 'yR, H is a Hubble parameter, and the dot means the derivative 
with respect to the time. The function m(r) [23] shows the deviation from 
the ACDM model. Using (21), (22) equations of motion (in the absence of 
the radiation, = 0 ) can be written in the form of autonomous equations 
[25] . The critical points for the system of equations can be investigated by the 
study of the function m{r). The plot of the function m{r) is given in Fig. 8 . 
The de Sitter point Pi [23] (in the absence of radiation, 0:4 = 0) corresponds 
to the parameters xi = 0 , X 2 = — 1 , X 3 = 2 (P = 0 ,r = — 2 ). One can check 
with the aid of Eqs.(6),(21) that Pi corresponds to the constant curvature 
solutions. The parameter of matter energy fraction, Om, and the effective 
equation of state (EoS) parameter, Wgff, are 


Hm = 1 - xi - 0)2 - X 3 = 0, tCgg = -1 - 2H/ (3P^) = -1. (23) 


These parameters correspond to DE. The constant curvature solution x ~ 
0.965 (a = 0.9) gives unstable the de Sitter space because 1 < m{r = —2) 


For the critical point P 5 ( 0:3 = 1/2), m ~ 0, r —1, and EoS of a matter 
era is = 0 (a = Then we have a viable matter dominated epoch 
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Figure 8 : The function m{r) 


prior to late-time acceleration [23]. The equation m = —r —1 has the solution 
m = 0 , r = — 1 , R = 0, corresponding to the point P 5 . One can verify with 
the help of Eq.(22) that m'{r = —1) = 0 (see also Fig. 8 ). Therefore the 
condition m'(r = — 1 ) > —1 is satisfied and we have the standard matter era 
As a result, there is the correct description of the standard matter era 
in the model under consideration. 


6 Conclusion 

The new model of modihed F{R) gravity is suggested and investigated. This 
is the effective gravity model which can describe the evolution of universe. 
The model gives the detailed account of the acceleration of universe corre¬ 
sponding to the de Sitter space and constant curvature solutions. It was 
shown that the de Sitter space-time is unstable because it corresponds to 
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the maximum of the effective potential. We have considered both the Jordan 
and Einstein frames and obtained the potential and the mass of the scalaron. 
It was demonstrated that flat space-time is stable. The slow-roll parameters 
e, r] and the e-fold number of the model were calculated. We demonstrate 
that the model can give the reasonable age of the inflation. The analysis of 
critical points of autonomous equations shows that the standard matter era 
exists and the necessary conditions for the standard matter era are satished 
in the model considered. Thus, the model can be alternative to GR, and can 
describe early-time inflation. 
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